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Abstract 

Aging in complex systems is studied via the sandpile model. Relaxation 
of avalanches in sandpiles is observed to depend on the time elapsed since the 
beginning of the relaxation. Levy behavior is observed in the distribution of 
characteristic times. In this way, aging and self-organized criticality appear to 
be closely related. 
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I. Introduction 



In the current literature of complex systems, aging is frequently reported 
as a phenomenon observed in. ajoki^lasses and other disordered systems, rather 
than in biological systemsj!lu'cK3 33 though some works have been devoted to 
biological aging of populations. uU Relevant contributions have been brought 
presenting this phenomenon from a general perspective, revealing that tempo- 
ral autocorrelation functionsJn the Bak-Sneppen model display aging behavior 
similar to glassy systemsHEij To our knowledge those papers give, up to now, 
the more relevant contribution to the relation between SOC and aging. 

Up to our knowledge, a detailed study of individual aging from the viewpoint 
of self-organized critical systems has not yet been performed. 

It is well known that biological aging manifests itself in an individual as 
the slowing down of many processes as, e.g., slower growing of tumors and 
manifestation of senescence in the slowing down of reflex behavior. Due to its 
characteristics, relaxation phenomena in complex systems are good candidates 
to the study of aging. 

About 150 years ago, Weber and Gauss carried out a simple experiment 
demonstrating that relaxation in complex systems is not exponential. Investi- 
gating the contraction of a silk thread they found that it does not contract so 
quickly but it relaxes slowly following a power law t~ a . This behavior is not 
particular of mechanical relaxation, it has been also observed in experiments of 
magnetic relaxation in spin glasses and high critical temperature superconduc- 
tors, transient curreni_.measurements in amorphous semiconductors, dielectric 
relaxation, and more.M 

In relaxation dynamics, aging means that the properties of a system depends 
on its ageJia For instance, consider a glass quenched at time t = below its 
glass transition temperature under an external stress. At time t = t w the stress 
is released. If the system were near equilibrium, then its response measured at 
certain time t > t w will be a function of the difference t — t w and, therefore, 
normalized responses taken at different initial times t w will collapse in a single 
curve. However, this is not the behavior in glass materials, and, as was pointed 
above, in live beings. 

Aging is a consequence of the nonequilibrium dynamics and may be consid- 
ered as a characteristic of the dynamics of complex systems far from equilibrium. 

In the following sections we demonstrate that these manifestations of indi- 
vidual aging are present in sandpile avalanches modeled- with a Bethe lattice, 
particularly, the process of relaxation. Contrary to,oli9 we do not consider 
deterministic "toppling" (like that present in the Bak-Sneppen model) but a 
probabilistic one, and avalanches are not regarded as infinite, as is,EI but the 
finite size of the system (sandpile) is explicitly included. 
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II. Sandpile model 



Sandpiles seem to be simplest systems which lead to complex behavior and 
are a paradigm for the study of all phenomena manifested in complex systems 
out of equilibrium, like relaxation under external perturbations. They have been 
taken as the paradigm of self-ocganization since the introduction of these ideas 
by Bak, Tang, and Wiesenfeld.ll3 The evolution of a sandpile, under an external 
perturbation has been extensively studied experimentally ,EJll3 theoretically ,ll3EZl 
and by computer simulations .La In general, those works analyze the sandpile 
dynamics close to the critical state, i.e., close to the critical angle. 

If a sandpile is subjected to an external perturbation, like, e.g., small am- 
plitude vibrations, large amplitude grain motions are rare events but from time 
to time, a grain may jump from its quasiequilibrium position. If a surface grain 
jumps, it will fall through the slope of the pile until it collides with another 
sand grains down the slope. After this collision the initial grain may be trapped 
with those grains or some of them may fall through the slope. If the last pos- 
sibility happens then each grain "surviving" the collision will fall through the 
slope as the initial grain did. This image of an avalanche as an initial object 
that consecutively drags another resembles a br-anching process for which the 
Bethe lattice representation seems to be natural.tDLa The representation of the 
avalanche dynamics as a Bethe lattice is a mean field approach to the problem. 
It neglects correlations between branches. Notwithstanding, if the pile angle 
is below is critical value the avalanches will be rare events and, in case of oc- 
currence, will be very sparse. Thus, such an approximation will be acceptable 
for analyzing the long time relaxation of the pile angle below its critical value, 
which is the subject of this paper. Away from the critical point the mean field 
approximation works quite well. 

Let us represent the avalanche as a cascade in the Bethe lattice as follows. 
Firstly, we start with a single node, which could represent in this case a grain. 
In a further step F will emerge with probability p(9), depending on the pile 
angle 9. This operation of generation of F identical particles starting from one 
is repeated in the next step to each node of the new group, and so on. If the 
percolation process overcomes a given length (that of the border of the pile) 
those nodes beyond the limit constitute the avalanche. If it does not, there 
would be no avalanche since the cascade was stopped before reaching the base 
of the pile (frustrated avalanche). By avalanche size we take the number of 
nodes, of the corresponding Bethe lattice, in the last step. 

After an avalanche the sandpile autoorganizes itself with the new number of 
grains (i.e., a new slope is calculated with the remaining grains). The occurrence 
of avalanches will carry as consequence a decrease in the number of grains in 
the pile and, therefore, a decrease of p(0). Each time an avalanche occurs, 
the occurrence of a new avalanche is less probable. Thus, to characterize this 
feedback mechanism a relation between p and the number of grains in the pile 
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N is needed. 

The drag probability p(9) is a function of the slope. Its value is determined 
by the competition of two contrary forces: the gravity, which conspires against 
the stability of the slope, and the friction which favors the slope stability. Since 
the slope forms an angle 9 with the horizontal plane the component of the 
gravity force in the slope direction will be larger with the increase of this angle, 
varying from zero to a maximum value when 9 goes from zero to tt/2. Therefore 
it is plausible to assume that the contribution of the gravity and, therefore, the 
tendency of falling down the slope is proportional to sin#. On the contrary, the 
resistance to this tendency given by the static friction decreases with decreasing 
the pile angle, varying from a maximum value to zero when 9 goes from zero 
to tt/2. Hence, it is also plausible to assume that the resistance to the falling 
down is proportional to cos 9. 

Thus the slope dependence of p can be expressed through the ratio of both 
tendencies sin 9/ cos 9 = tan 9. Based on this hypothesis we may propose the 
exponential relation 

p{9) = exp(-A/tan0) , (1) 

where A is a parameter determined by the gravitational field, the friction, and 
vibration intensity. Notice that p(Q) = and p{ir/2) = 1. This selection for 
p{9) can not be regarded as a sophisticated trick or something imposed ad hoc 
to the model in order to obtain the desired results, since it can be verified that 
any function satisfying the very general conditions already stated is adequate 
to our model. Incidentally, it must be said that this reveals the robustness of 
the model. 

On the other hand, the number of grains in a pile with slope angle 9 is given 

by 

N = Bx 3 tan 8, (2) 

where B is a geometrical factor and x = D/d is the ration between a character- 
istic size of the pile base D and sand grain d. Combining equations (|l|) and (||) 
it is obtained 

p(N) = exp(-ca; 3 /^V) . (3) 

where c = AB. Thus Equation [l] relates the dragging probability with the 
number of grains in the pile. 

The parameter c should not depend on the size of the system. It must be a 
function of parameters describing the vibrations, amplitude a and frequency to, 
and of the gravitational field acceleration g. The only nondimensional combina- 
tion of these magnitudes is given by the ratio between the vibration acceleration 
aoj 2 and g and, therefore, c = c(aio 2 /g). This conclusion, obtained from dimciL. 
sional analysis, is corroborated by experiments on sandpiles under vibrational! 
which show that the ratio auj 2 /g is the relevant parameter. 
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For a given value of c and x , there is a critical number of grains in the pile 
N c . This value can be found recalling that in the Bcthc lattice the critical value 
of p for percolation exists is 1/F, resulting 

N c = cx 3 /\nF . (4) 

In a static sandpile (i.e., no vibrations) thiji-value corresponds with the number 
of grains in the pile at the angle of reposeO 

Another magnitude of interest is the penetration length, the number of steps 
n(6) in the Bethe lattice for an avalanche to take place. This magnitude must 
be proportional to the length of the pile slope and, therefore, should be given 
by the expression 

n(6)=x/cos6 (5) 

Notice that the possible existence of a geometrical factor in equation || may be 
absorbed in x, redefining B in Equation ^. Moreover, the relation between n 
and N may be easily obtained using Equation |^. 

Equations || and [| relate the parameters of the pile with those of its Bethe 
lattice representation. -They were obtained here in a different way than that pro- 
posed by Zaperi et a/.Lj Other dependencies between the dragging probability 
and the number of grains in the sandpile may be proposed. Notwithstanding, 
as it is discussed below, the precise form of this functional dependence is not 
relevant. 



III. Simulations and results 



The numerical experiment of relaxation is performed as follows. We start 
with a certain number of N. In a first step, we test if an avalanche takes place 
using the Bcthc lattice representation. If it does, then we simulate the process 
and then recalculate the value of p(N) by simply substracting the size of the 
avalanche to N and using Equation g. The size of the avalanche is the number 
of nodes generated in the Bethe lattice which surpass the size of the pile whose 
measure is given by (|^). Then this step is repeated again and again, whereas 
avalanche sizes and times are registered. 

If avalanches are considered as instantaneous the number of steps is a mea- 
sure of time. This approximation is valid for low vibration intensities. In this 
case grain jumps which trigger avalanches are rare events and, therefore, the 
time between two successive grain jumps will be much larger than the duration 
of avalanches. 

To describe the behavior of our system, we define the relaxation function 
4>{t + t w ,t w ) of the sandpile as: 

, . , s N(t + t w ) 
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Figure 1: Aging. Normalized relaxation function <f)(t + t w ) = N(t + t w )/N(t w ) 
taken at three different initial times t w . The long time tail exhibited in the 
figure has an exponent a ~ 0.089. 

where we include the time dependence of the number of grains in the sand- 
pile. t w is the "waiting time", i.e., the instant at which we begin the coun ting 
of the number of grains since the start of the relaxation, (avalanche), as in.BEJ 

The existence of aging in our model is illustrated in Figure 1. We have 
plotted the normalized relaxation <j>{t + t w ,t w ) at different ages (steps) of the 
system (simulation) by taking as initial time t — t w for x = 10 and c = 1. 
As it can be seen the relaxation becomes slower with increasing the age. of the 
system, in agreement with experiments in structural and spin glassea!3 (and, 
unfortunately, with human life). The system exhibits a delay in the relaxation, 
after which the relaxation function decreases as a power law (long time tail) with 
an exponent a ~ 0.089. Thus, since the angle relaxation becomes slower with 
the age of the pile then it will never reach an equilibrium angle and, therefore, 
properties like translational invariance and the fluctuation dissipation theorem 
do not hold £3 

Associated with these slow relaxation dynamics and aging phenomena we 
expect to observe a wide distribution of time between avalanches. As time in- 
creases the time between two consecutive avalanches At becomes larger, because 
after an avalanche the occurrence of a new avalanche becomes smaller. There- 
fore, it is expected that the mean time between avalanches diverges as t — > oo. 
Hence, the distribution of time between avalanches n(At) should satisfy the 
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Figure 2: Distribution of time between avalanches. The dashed line is a power 
tail with exponent — 1 — a , a being the exponent of the long time tail shown 
in Figure 1. 

asymptotic behavior for large At 

ra(At) ~ Ar 1 -?, (7) 

with < j3 < 1. 

This hypothesis is confirmed in our simulations. Figure 2 shows the distri- 
bution of time between avalanches for x = 10 and c = 1. It is approximately 
constant for small values and then it decays following a power law in more than 
two decades. The plateau at small values of At is associated to the rapid decay 
observed at short times while the tail for large At should be related to the long 
time tail. Here we observe that (3 ~ a ~ 0.089 Thus, there should be some 
connection between the distribution of time between avalanches and the long 
time relaxation. 

Finally we want to mention that these simulations were also carried out 
assuming other functional dependences between the dragging probability and 
the number of grains in the pile p( N) . In all cases the results where qualitatively 
similar to those presented here using equation [| reflecting that the precise 
dependence is not relevant. 
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IV. Conclusions 



A sandpile model for aging was presented revealing that this phenomenon 
is manifested in relaxation of sandpile avalanches. Biological systems and indi- 
viduals show similar behavior. Phenomena related to some kind of relaxation 
in live beings deserve more quantitative research. The introduction of a Bethe 
lattice representation for the avalanches and a feedback mechanism describes 
quite well the principal features of the relaxation in sandpiles under low inten- 
sity vibrations. The proposed representation leads to long time tails relaxation, 
aging and Levy (fractal) distributions of time constants, which are characteristic 
properties of the dynamics of complex systems out of equilibrium. 
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